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ABSTRACT 

We report the detection of an extremely bright (~37 mJy at 1100 /xm and ^^91 mJy 
at 880 Aim) submilhmeter galaxy (SMG), AzTEC-ASTE-SXDFllOO.OOl (hereafter 
referred to as SXDFllOO.OOl or Orochi), discovered in 1100 /im observations of the 
Subaru/XMM-Newton Deep Field using AzTEC on ASTE. Subsequent CARMA 1300 
/zm and SMA 880 /im observations successfully pinpoint the location of Orochi and 
suggest that it has two components, one extended (FWHM of ^ 4") and one compact 
(unresolved). Z-Spec on CSO has also been used to obtain a wide band spectrum 
from 190 to 308 GHz, although no significant emission/absorption lines are found. 
The derived upper limit to the line-to-continuum fiux ratio is 0.1-0.3 (2 a) across the 
Z-Spec band. 

Based on the analysis of the derived spectral energy distribution from optical 
to radio wavelengths of possible counterparts near the SMA/CARMA peak position, 
we suggest that Orochi is a lensed, optically dark SMG lying at z ^ 3.4 behind a 
foreground, optically visible (but red) galaxy at z 1.4. The deduced apparent (i.e., 
no correction for magnification) infrared luminosity (im) and star formation rate 
(SFR) are 6 x 10-^^ and 11000 Mq yr^^, respectively, assuming that the Lir is 
dominated by star formation. These values suggest that Orochi will consume its gas 
reservoir within a short time scale (3 x 10^ yr), which is indeed comparable to those 
in extreme starbursts like the centres of local ULIRGs. 

Key words: galaxies: high-redshift - galaxies: starburst - galaxies: ISM - submil- 
limeter 
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1 INTRODUCTION 

Recent rapid advancements in wide and deep surveys at mil- 
limeter/submillimeter (mm/submm) wavelengths have led 
to successive discoveries of numerous mm/submm-bright 
galaxies (SMGs) in the early universe. 

After the pioneering works performed wit h the SCUBA 
l|Cunningham et al ]|l994l:lHolland et al.l | l999l) on the JCMT 



15 m telescope (e . g.TISmail e t al.' '19971: iHuehes et all 1 19981 : 
lEales et"al] 1 19991 : llvison et a l. 2000), numerous mapping 
surveys at A = 850 /jm have been conducted toward blank 
fields (e.g.. lBorvs et al.ll2003l : rCoppin et al .ll2006l:IPop e et all 
2006h and over-den se regions (e.g.. [Stevens et al.l 



20031 : 



Knudsen et "aLll2006l ). MAMBO jKrevsa et alj|l998ll on the 



IRAM 30 m telescope and BOLOCAM ijLaurent et al.ll2005l ) 
on the CSO 10 m have also been used to obtain 1200 
/xm/llOO nm image s of blank fields (Gre ve et al. 2004, 2008; 
iBertoldi et al.ll2007l: iLaurent et al.ll2005l l and over-dense re- 
gions (e.g. jGreve et al.ll200'i1 )" 

A ne w mm-wavelength b olometer camera, AzTEC for 
1100 ^m l| Wilson et al.ll2008al l. was mounted on JCMT and 
produced wide-area (~ a few 100-1000 arcmin' ^ scale) im- 
ages of well studied fields such a s COS MOS (jScott et al.l 
120081 ) ■ GOODS-N llPerera et al.l |2008| ). Lockman hole, 
SXDF llAustermann et all |2010| ). and over- dense regions 
ijWardlow et al.ll2010h . 

New submm wave telescopes in northern Chile, i.e., 
ASTE 10 m and APEX 12 m telescopes, are now also 
equipped with bolometer camer as, i.e., AzTEC fo r 1100 
/xm and LABOCA for 870 fj.m jSiringo et aLlbOOSl ). Ow- 
ing to the very suitable atmospheric conditions of the site, 
these telescopes routinely ob tain wide-area images of var- 
ious fields su ch as SSA22 jTamura et all [2009.) . ECDF- 
S/GOODS-S llWeifi et al ]|2009l: IScot t et al.l l2010l). A DF-S 
(jHatsukade et al. 12011 ), a proto-cl uster ( Beelen et al.ll2008l ) 
and a cluster ( Wilson et al.ll200"8bh . A wider area (~10 deg^ 
scale) short su bmm survey has been conducted with BLAST 
(jPevlin et al.i i2009). and ~100 deg^ scale rnm/su bmm sur- 
veys a re now coming using SPT JVieira et al.|l2010l ') and Her- 
schel (jAmblard eralll2010l : lEales et al.ll20ld '). 

One of the important findings from these recent 
mm/submm surveys is the detection of ultra-bright pop- 
ulations of SMGs and it is predicted that this ultra- 
bright population is probably l ensed by foreground clus- 
ters a nd /or massive galaxies (e.g.. lBlciinll 19961 : iNegrello et al.l 
I2OO7I ). For instance, an ultra-bright mm/submm galaxy, 
MM J065 837-5557.0, with a fl ux density of ~ 20 mjy at 
1100 um dWilson et al] 12008^ ) and ~ 48 mJy at 870 ^m 
(jjohansson et al.ll2010l ) has been detected near the center of 
the Bullet cluster {z = 0.297). This source, whose bright- 
ness appears to be highly boosted (> 20-75) by gravita- 
tional lensing, is a luminous infrar ed galaxy (LIRG) be- 
hind the Bullet clust er at z ~ 2.79 (jjohansson et al]|2010l : 
IConzalez et al.l I2OI0I ). A similar but brighter source has 
also been reported toward the cluster MACS J2135-010217 
{z = 0.325). SMM J2135-0102 has an 870 fim flux density of 
106 mJy, and is an ultra-luminous IR galaxy (ULIRG) at a 
spectroscopically confirmed red shift of 2.33 and with an am- 
plifica tion factor of ~33 ( Swinb ank et al.|[201(]| : llvison et al.l 
l2010al ). Owing to the strong magnification of the gravita- 
tional lens, these lensed SMGs provide a unique opportu- 
nity to understand the physical properties of extreme star- 



formati on in the early universe even with existing tele- 
scopes (jSwinbank et all 120101 ). Furthermore, SPT surveys 
have shown that ultra-bright submm/mm galaxies exist in 
some surface density; 20 dust-dominated SM Gs above 10 
mJy at 1.4 mm were detected in 87 deg^ area (jVieira et al.l 
|201(]| ). The 14.4 deg survey with Herschel, as a part of the 
H-ATLAS project, uncovered 11 bright 500 /^m sources (> 
100 mJy at 500 fim) within the survey area, and 5 of 11 
have been identified as lensed, dusty starburst galaxies at 
z = 1.6 — 3.0 (NcgrcUo et al. 2010). Large mm/submm sur- 
veys enable us to find these new population. 

Here, we report a serendipitous detection of an ultra- 
bright SMG in the Subaru/XMM-Newton Deep Field 
(SXDF), termed as AzTEC- ASTE-SXDFllOO.OOl (hear- 
after referred to as SXDFllOO.OOl) or Orochfl during a 
course of wide and deep 1100 /im surveys using AzTEC 
mounted on ASTE. 

This paper is organized as follows. The AzTEC on 
ASTE detection is reported in section 2, and subsequent 
CARMA and SMA observations are described in section 
3, along with the multi-wavelengths images from opti- 
cal/infrared to radio. Spectroscopic observations using Z- 
Spec on CSO are discussed in section 4. The modeling of 
the SED and derived physical properties, including discus- 
sions on the source size/structure and implications for star 
formation properties of Orochi are given in section 5. 

Throughout this paper, we adopt a cosmology with den- 
sity parameters Qa = 0.7 and Q.m = 0.3 and the Hubble 
constant Ho = 70 km s~^ Mpc~^. 



2 AZTEC/ ASTE 1100 fiM OBSERVATIONS 
AND RESULTS 

2.1 AzTEC/ ASTE observations 

We conducted 1100 /im imaging observations o f a cluster of 
LyQ Emitters (LAEs) at z 



5.7 (Clump B; Ouchi et al 



12005") in SXDF using the AzTEC camera (jWilson et aL 
2008a) mounted on ASTE (Ezawa et al. 2004, 2008), from 
November 26 to December 21, 2008. The observations 
were carried out remotely from the ASTE operation rooms 
through the network observation system N-COSMOS3 de- 
veloped by the National Astronomical Observatory of Japan 
(NAOJ) (jKamazaki et al. 2005). The full width at half max- 
imum (FWHM) of the AzTEC beam on ASTE is 30" at 1100 
^m, and the field of view of the array is roughly circular with 
a diameter of 8'. During the ASTE 2008 observation run, 117 
of 144 AzTEC detectors were operational. 

We imaged a ~6.6' diameter circular field of Clump 
B, centered at R.A. (J2000) = 02*^ 18" 19.60", Dec. 
(J2000) = -5°32'52 ".0. We used a Lissajous scan pattern 
( Wilson et alll2008bl ) in order to maximize the observation 
efficiency. We selected a maximum velocity of 300" s~^ in or- 
der to mitigate low-frequency atmospheric fiuctuations. We 
obtained a total of 39 individual observations for Clump 
B, taking ~40 min for each observation. The atmospheric 
zenith opacity at 220 GHz was r220GHz = 0.017-0.100 as 
monitored with a radiometer at the ASTE telescope site. 



^ a Japanese word referring to an ancient Japanese legendary 
monster. 
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Table 1. AzTEC on ASTE observations. 



Parameters 


Values 


Observation date 


Nov. 26 - Dec. 21 2008 


^Vave 1 e n g t h. / Fr e Q u e ncy 


1100 /im/270 GHz 


j_) diiti w m Lii 






^AA. /'fntfin 117 f^nnor::ifinnan 


Beam size (FWHM) 


30" 


Field center (J2000) 


R.A. = 02'' 18™ 19.59" 




Dec. = -05°32'52".00 


Field area 


136 arcmin^ (50 %-coverage region) 


Map noise level 


0.6-1.0 mjy 


Pointing source 


J0132-169 


Flux calibrator source 


Uranus, Neptune 


Opacity at 220 GHz 


0.017-0.100 



Uranus or Neptune were observed at least once a night 
in order to measure each detector's point spread function 
(PSF) and relative position and to de termine the flux con - 
version factor for absolute calibration (| Wilson et al.l[2008al ). 
Pointing observations with the quasar J0132-169 were per- 
formed every 2h across observations for Clu mp B; the re- 
sultan t pointing accuracy is better than 3" (| Wilson et al.l 
l2008bl ). A pointing model is devised by interpolating these 
pointing data temporally and is applied to the astrome- 
try for correction of pointings. Observational information 
is summarized in Table [T] 



2.2 AzTEC Data Reduction 

The data were reduced using the AzTEC data reduction 
pipeline written in the Intera ctive Data Langua ge (IDL) in 
a manner similar to that in IScott et al.l (|2008l ). Here, we 
provide a brief summary of the process and point out the 
difference. The data were divided into 15 s intervals of time- 
series data. Spikes were then removed from the time-series 
data. A principal component analysis (PCA) method was 
used to subtract the sky emission, and the effect of PCA 
cleaning on the point source response is traced by reducing 
synthetic, noiseless time-series data with a simulated point 
source using the same approach that employed for the ac- 
tual data. This 'point source kernel' is indicative of the effect 
of PCA on the point source profile and flux attenuation of 
real point sources in the map and we correc ted flux densi- 
ties o f AzTEC sources based on this kernel (jPownes et al.l 
[20l3). The cleaned time-series data were projected into a 
map space using 3" x 3" pixels, and the 39 individual ob- 
servations were co-added into a single map by weighted av- 
eraging. We also created 100 synthesised noise realizations 
by randomly multiplying each 15 s time-series interval by 
± 1 (similar to t he 'sc an-by-scan' jackkniflng technique used 
in I Scott et al.l (|2008l )). These noise realizations are free of 
astronomical signals, including the signals from confused 
sources in the map. The co-added map and the 100 noise 
realizations were optimally flltered for the detection of point 
sources. 




Right Ascension (J2000) 



Figure 1. The AzTEC/ASTE 1100 /xm image of Orochi. The 
contour levels are 5, 10, 15, 20, 25, 30, 35, 40, 45, 50 and 55 a, 
where 1 a = 0.65 mJy beam~^. Negatives are shown in white 
contours. 

2.3 AzTEC 1100 results 

The achieved average noise level of the resulting map is 
0.60-1.00 mJy over a 136 arcmin^ area. We found an ul- 
tra bright source, Orochi, at a(J2000) = 02''18™30.68'' and 
5(J2000) = -05°31'31".37. The flux density at 1100 fim is 
37.28 ± 0.65 mJy. The 1100 fim image of Orochi is shown 
in Fig. 1. The source size is consistent with the point source 
kernel or PSF of the AzTEC/ASTE observations. Hence, 
the point source kernel is understood as the AzTEC/ASTE 
beam. The estimated position errors from AzTEC signal- 
to-noise ratio and the beam size are 0".9 (1 a) and l".4 (2 
a). 



3 INTERFEROMETRIC IDENTIFICATION 
AND MULTI WAVELENGTH PROPERTIES 

We performed interferometric imaging with CARMA and 
SMA of Orochi to constrain its position. With the reflned po- 
sition estimate, we searched the literature and the archives 
for multi-wavelength data. 

3.1 CARMA 1300 observations and results 

3.1.1 CARMA 1300 fim Observations & Reduction 

We conducted follow-up observations of Orochi on August 
16, 2009, using the D configuration with 15 antennas of 
the CARMA. The phase centre was a(J2000) = 02''18™21" 
and (5(J2000) = -05°3l'3l". The projected baseline lengths 
ranged from 10 m to 108 m. The CARMA correlator was 
configured to cover a 1.5 GHz width in each sideband, yield- 
ing a total bandwidth of 3 GHz for continuum observations 
after adding the two sidebands. The centre frequency of the 
receivers was tuned to 230 GHz. We observed a bright QSO, 
J0108-I-015 (18.8° away from Orochi) as a visibility calibra- 
tor and Uranus as a flux calibrator. To obtain an empirical 
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Table 2. CARMA observations. 



Parameters 


Values 


Observation date 


AiiD- 16 2009 


Wavelengtli/Frecjuency 


1300 /^m/232 GHz 


Bandwidth 


1.5 GHz 


Phase center (J2000) 


R.A. = 02''18"'21'' 




Dec — — 05°3T31" 


Phase cahbrator 


J0108+015 


Flux calibrator 


Uranus 


Array configuration 


D configuration 


projected baseline 


10-108 m 


Primary beam 


35". 4 (FWHM) 


Synthesized Beam size 


3". 2 X 2".l {P.A.= -16.4°) 


Map noise level 


1.3 mjy 


Opacity at 230 GHz 


0.1-0.3 



upper limit on the systematic position error induced by base- 
line errors, we observed a radio galaxy, J0241-082 (25° away 
from the visibility calibrator; the distance is approximately 
1.33 times that between Orochi and the visibility calibra- 
tor) during a track. We observed a bright QSO, 3C84, as a 
bandpass calibrator. The raw CARMA data wer e calibrated 
and i maged with natural weight using MIRIAD jSault et al.l 
Il995f ) . We found that the quality of the data from the upper 
side band (USB) was significantly poorer than that of the 
lower side band (LSB) data. Therefore, we used only the 
LSB data for pinpointing the position of Orochi and subse- 
quent analysis. Observational information is summarized in 
Table 2. 



3.1.2 CARMA 1300 results 

We find a source with 11 a significance at a(J2000) — 
02''18™30.67'' and (5(J2000) = -05°3l'3l".42 (Fig. [2}. The 
CARMA 1300 pm source position coincides well with that 
of the centroid of the AzTEC/ASTE 1100 fim source. The 
derived source position is shown in Table [3] 

To evaluate the reliability of the astrometry in our 
CARMA observations, we estimate the statistical errors due 
to the noise obtained in fitting a point source to the cali- 
brated visibilities and systematic errors due to uncertain- 
ties in i nterferometer baselin es. Details of these analysis 
followed lYounger et all (|2007| ). The derived statistical er- 
rors are shown in Table |31 The systematic astrometry error 
caused by baseline length uncertainties is estimated as ~0".5 
based on the CARMA image of a known radio source J0241- 
082. The estimated statistical positional errors for RA and 
Dec are 0".09 and 0".12. 

3.2 SMA 880 /im observations and results 

3.2.1 SMA 880 fim Observations & Reduction 

We conducted follow-up observations of Orochi on December 
10, 2009, using the compact configuration with 8 antennas 
of the SMA. The phase center was q(J2000) = 02''18™22^ 
and (5(J2000) = -05°31'35". The projected baseline lengths 
ranged from 7 m to 70 m. The SMA correlator was equipped 
with 4.0 GHz in each sideband, providing a total of 8 




o 02''18'"31f5 31?0 30?5 30f0 29?5 



Right Ascension (J2000) 

Figure 2. The CARMA 1300 fim image of Orochi. The contour 
levels are 3, 6, and 9 (t (1 ct = 1.3 mJy). Negatives are shown 
in white contours. The synthesized beam is shown on the lower 
left. Note that no correction of the primary beam attenuation was 
made in this map. The larege solid yellow circle shows the beam 
size of AzTEC/ASTE (30"). The small solid white circles show 
the estimated position errors; the inner one shows 1 a error (0".9) 
and the outer one shows 2 a error (1".4). 



Table 4. SMA observations. 



Parameters 


Values 


Observation date 


Dec. 10, 2009 


Wavelength /Frequency 


880 ^lm/ 340 GHz 


Bandwidth 


8 GHz 


Phase center (J2000) 


R.A. = 02''18'"22'' 




Dec. = -05°3T35" 


Phase calibrator 


J0132-169 


Flux calibrator 


Uranus 


Array configuration 


compact configuration 


projected baseline 


7-70 m 


Primary beam 


32". 4 (FWHM) 


Synthesized Beam size 


3".0 X 1".9 (P.A.= 8.1°) 


Map noise level 


2.2 mJy 


Opacity at 225 GHz 


0.1-0.3 



GIfz bandwidth for continuum observations. The centre fre- 
quency of the receivers was tuned to 340 GHz. We observed 
a bright QSO, J0132-169 (16.0° away from Orochi) as a vis- 
ibility calibrator and Uranus as a flux calibrator. To obtain 
an empirical upper limit on the systematic position error in- 
duced by baseline errors, we observed a radio galaxy, J0238- 
166 (22.7° away from the visibility calibrator; the distance is 
approximately 1.4 times that between Orochi and the visi- 
bility calibrator) during a track. We observed a bright QSO, 
3C454.3, and Callisto as bandpass calibrators during the 
track. 

The raw SMA data we re calibrated using the MIR 
package (|Scoville et al.l Il993l ) . Imaging was carried out in 
MIRIAD with natural weighting. Observational information 
is summarized in Table 4. 
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Table 3. Peak positions and their errors for VLA 20 cm, CARMA 1300 fim, AzTEC 1100 fim, SMA 880 fim, and Subaru z'-band images 
of Orochi. The measured offsets from the z'-band peak position are also listed. 











a{a) 


a{5) 




(t'(5)'' 


z'-band offset cr{a) 


z'-band offset cr(S) 


image 


a(J2000) 


(5(J2000) 




(arcsec) 


(arcsec) 


(arcsec) 


(arcsec) 


(arcsec) 


(arcsec) 


VLA 


02'»18™30.67'' 


-05°31'31" 


48 


0.07 


0.09 


0.26 


0.23 


-0.18 


0.32 


CARMA 


02'^18™30.67'' 


-05°31'31" 


28 


0.09 


0.12 


0.27 


0.24 


-0.18 


0.52 


AzTEC 


02''18™30.67= 


-05°31'30" 


97 


0.9 


0.9 


0.9 


0.9 


-0.18 


0.83 


SMA 


02'"18™30.68'' 


-05°31'31" 


68 


0.04 


0.06 


0.25 


0.22 


0.00 


0.12 


z'-band 


02'"18'"30.68'' 


-05°31'31" 


80 


0.25 


0.21 











" denotes the combined (z' plus radio or submm/mm) 1 cr uncertainty in position. 
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Figure 3. The SMA 880 ^m image of Orochi. The contour levels 
are -3, 3, 6, 9, 12, 15, 18, and 21 tr (1 tr = 2.2 mjy). Negatives are 
shown in white contours. The synthesized beam is shown on the 
lower left. Note that no correction of the primary beam attenua- 
tion was made in this map. The large solid yellow circle shows the 
beam size of AzTEC/ASTE (30"). The small soUd white circles 
show the estimated position errors; the inner one shows 1 cr error 
(0".9) and the outer one shows 2 cr error (1".4). 



3.3 Multi-wavelength Data 

We find a multi-wavelength counterpart at the 
CARMA/SMA position using optical, near-infrared 
(NIR), mid-infrared (MIR), and radio wavelengths from 
archival data. The derived flux and corresponding images 
are displayed in Table 5 and Fig. 4, respectively. Here, 
we briefly describe the archival data. Optical data in 5 
broad-bands, B, V, Rc, z' , and i'- band, were provided b y 
the SXDS project data release 1 (jPurusawa et al.l |2008| ) . 
A iifl-band image was obtained from the SMOKA archive 
in Japan, if -band and J-band imag es were provided b y 
UKIDSS/UDS surveys data release 1 (jWarren et al.ll2007l ). 
Infrared data in all four Spitzer IRAC bands (3.6, 4.5, 5.8, 
and 8 /im) are available from the SWIRE survey, and the 
MIPS band (24 ^m) is available from the SpUDS survey. 
These data were obtained from the Spitzer archive. Orochi 
was detected in a 20 cm VLA image a nd catalogued as VLA 
J021830-05315 (|Simpson et all |2006| ). An updated VLA 
image, taken in the A-configuration (Arumugam et al., in 
prep.), has been incorporated into the multi- wavelengths 
dataset. Unfortunately, the counterpart of Orochi is not 
isolated on the 24 /xm image . Other infrared bright source 
is located in south west of Orochi and it looks like an 
extended source in the optical wavelength images (Fig. |4]). 
Therefore, we regard the aperture photometry data in the 
SWIRE catalogue as an upper limit. The photometry data 
at 1000, 1100, 1200, 1300, 1400, and 1500 fim are obtained 
from Z-Spec observations. See section 4 for details. 



3.2.2 SMA 880 results 

We find a source with 21 a significance at the position of 
q(J2000) = 02''18'"30.67'' and 5(J2000) = -05°31'31".68 
(Fig. |3]). The SMA 880 fim source position coincides with 
that of the centroid of the AzTEC/ASTE 1100 nm source 
and CARMA 1300 fim source. The derived source positions 
of AzTEC/ASTE, CARMA, and SMA are shown in Table 
[3] along with the VLA source position. 

The systematic astrometry error of the SMA source 
caused by baseline length uncertainties was estimated as 
~0".l based on the SMA image of J0238-166. The estimated 
statistical positional errors for RA and Dec are 0".04 and 
0".06. 



4 Z-SPEC 1500 TO 1000 

ULTRA- WIDEBAND SPECTROSCOPY 

4.1 Observations and data reduction 

In order to determine the redshift of Orochi, we searched 
for redshifted molecular/atomic lines using Z-Spec, a single- 
beam grating spectrometer which disperses the 190-308 
GHz (or 15 70-970 nm) band acr oss a linear array of 160 
bolometers (jBradford et al.l |2009| . and references therein). 
The resolving power (R) runs from 250 to 300, from the 
highest end to the lowest end of the band. This gives a fre- 
quency resolution of ~ 700 MHz or a velocity resolution of 
~ 880 km s~^ at the centre of the band (~ 240 GHz). 

Orochi was observed for six nights during late Novem- 
ber to early December of 2009 using Z-Spec mounted on 
the Caltech Submillimeter Observatory (CSO) 10 m tele- 
scope. A traditional chop-and-nod mode, with the secondary 
chopping at 1.6 Hz and a nod period of 20 s, was adopted. 
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Figure 4. Multi-wavelength images of Orochi with SMA contours (3, 6, 9, 12, 15, 18, and 21 a). The size of each image is 10" X 10"; 
north is towards the top of the image and east is towards the left side of the image. From left to right and top to down: VLA 20 cm; 
CARMA 1300 Atm; SMA880 /xm; MIPS 24 ^tm; rgb image of IRAC/chl (blue), ch2 (green), chS and ch4 (red); rg image of WFCAM/J 
(green) and K (red); rgb image of SuprimeCam/i?, i' and z'; rgb image of SuprimeCam/M/j, B and V. 



Table 5. Flux densities of Orochi. 



Wavelength 


Camera /Telescope 


Magnitude/Flux 


Reference 


4.5 - 10 keV 


EPIC/XMM-Newton 


< 5 Xl0"^'5(erg/cm2/s)'" 


Ueda et al. ('2008) 


0.5 - 4.5 keV 


EPIC/XMM-Newton 


< 8 Xl0"i''(erg/cm2/s)<' 


Ueda et al. (2008) 


0.380 fim (ufl-band) SuprimeCam/Subaru 


< 26.76 {niAB)'' 




0.445 urn (B-band) 


SuprimeCam/Subaru 


25.86 ± 0.05 iniAB) 


Furusawa et al. (2008) 


0.551 Aim (y-band) 


SuprimeCam/Subaru 


25.51 ± 0.05 irriAB) 


Furusawa et al. (2008) 


0.659 fj,m (i?c-band) 


SuprimeCam/Subaru 


24.80 ± 0.05 (m-As) 


Furusawa et al. (2008) 


0.771 Atm (i'-band) 


SuprimeCam/Subaru 


23.68 ± 0.05 {niAB) 


Furusawa et al. (2008) 


0.922 ^lm (z'-band) 


SuprimeCam/Subaru 


22.88 ± 0.05 (m^s) 


Furusawa et al. ('2008) 


1.215 Aim (J-band) 


WFCAM/UKIRT 


0.0096 ± 0.00066 (mjy) 


Warren et al. (2007) 


2.179 Aim (ii'-band) 


WFCAM/UKIRT 


0.0224 ± 0.0017 (mJy) 


Warren et al. (2007) 


3.6 fim 


IRAC/Spitzer 


0.046 ± 0.0018 (mJy) 




4.5 /im 


IRAC/Spitzer 


0.057 ± 0.0017 (mJy) 




5.8 Aim 


IRAC/Spitzer 


0.075 ± 0.0088 (mJy) 




8 Aim 


IRAC/Spitzer 


0.110 ± 0.013 (mJy) 




24 Aim 


MlPS/Spitzer 


< 0.384 (mJy) 




70 Aim 


MlPS/Spitzer 


< 12 (mjy)'= 




160 Aim 


MIPS /Spitzer 


< 72 (mjy)'= 




880 A^m 


-/SMA 


90.7 ± 2.2 ± 20.7'* (mJy) 


This work 


1000 Atm 


Z-SPEC/CSO 


51.9 ± 1.2 (mJy) 


This work 


1100 Atm 


AzTEC/ASTE 


37.3 ± 0.7 (mJy) 


This work 


1100 Aim 


Z-SPEC/CSO 


38.4 ± 0.9 (mJy) 


This work 


1200 Atm 


Z-SPEC/CSO 


31.3 ± 0.8 (mJy) 


This work 


1300 Atm 


-/CARMA 


24.6 ± 1.3 ± 5.3"^ (mJy) 


This work 


1300 Atm 


Z-SPEC/CSO 


23.9 ± 1.0 (mJy) 


This work 


1400 Atm 


Z-SPEC/CSO 


19.5 ± 0.9 (mJy) 


This work 


1500 Atm 


Z-SPEC/CSO 


13.9 ± 1.2 (mJy) 


This work 


20 cm 


-/VLA 


0.238 ± 0.038 (mJy) 


Arumugam et al., in prep 


50 cm 


-/GMRT 


0.315 ± 0.120 (mJy) 


Dunne et al. (2009) 



" indicates the detection limit in Ueda et al. 2008 and corresponds to a confidence level of 
99.91%. 

indicates the detection limit with 2 cr and aperture size of 2" (diameter), 
indicates the upper limit with 3cr. 

indicates the uncertainty of uv model fitting for the extended structure. See section 5.1 for 
details. 
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Table 6. Z-Spec observations. 



Observation date Values 

Observation date Nov.-Dec, 2009 

Wavelength/Frequency 970-1570 fj.m/ 190-308 GHz 

Spectral resolution ~ 700 MHz at 240 GHz 

Flux calibrator Mars 

Averaged noise level ~ 2.7 mjy 

Opacity at 225 GHz 0.03-0.17 



The total integration time was 16 h, excluding bad scans 
(showing unreasonable baseline shape and/or offset). The 
absolute intensity scale was calibrated from observations of 
Mars, with a channel-to-channel correction based on quasar 
spectra measured during the observation run. The overall 
calibration uncertainties are estimated to be less than 10%, 
except at the lowest frequencies which are degraded by the 
wing of the 186 GHz atmospheric water line. The Z-Spec 
observations are summarized in Table 6. 



4.2 Z-Spec result 

The derived Z-Spec spectrum of Orochi is displayed in Fig. 
[S] Continuum emission is detected across the whole band, 
and the S/N ratio exceeds 10 above ~220 GHz. The averaged 
noise level of the spectrum is 2.7 mJy. We found a clear slope 
of continuum emission, which is attributed to the thermal 
dust emission. The best-fit continuum spectrum across the 
Z-Spec band is 

/ ly \3.00±0.11 

^^-(30.1 ±0.5) (^3^) mJy. (1) 

Considering absolute flux calibration errors of 10-20% in the 
mm wavelength band, the measured flux densities with SMA 
and CARMA coincide well with the Z-Spec continuum flux. 
We derived the average flux densities at 1000, 1100, 1200, 
1300, 1400, and 1500 /im by binning the Z-Spec spectrum 
into bands of width ~7 GHz. The results are listed in Table 

m 

We see no significant emission/absorption line feature 
across the band, and therefore we were unable to determine 
the redshift of Orochi from these Z-Spec data. For a more 
quantitative analysis, we use the redshift-finding algorithm 
described in Lupu et al. (in prepration.), which computes 
the average S/N ratio from ah CO, [CI], [CII], and [Nil] 
transitions falling into the Z-Spec bandpass as a function of 
redshift. The resulting redshift probability function is very 
noisy, and the average S/N ratio at all redshifts is < 3 sigma, 
demonstrating the lack of significant line emission in these 
data. 

The upper limits on the expected CO emission lines and 
their implications are discussed in section 5. 



5 DISCUSSION 

In this section, we discuss the properties of Orochi, an 
uhra- bright SMG identification in AzTEC/ASTE surveys 
and confirmed by subsequent CARMA/SMA observations. 



First, the source structure is; the discovery of the spatially 
extended submm/mm bright component accompanied with 
an unresolved compact bright source is discussed. Then, de- 
tailed modeling of SEDs as well as the derivations of physi- 
cal quantities are presented. We suggest that an ultra-bright 
SMG at a redshift of ~3.4 has been discovered just behind 
a red foreground galaxy at 2: ~ 1.4 and that Orochi is likely 
to be lensed by the foreground galaxy. 

The deduced star formation properties and internal 
structures of Orochi are also discussed, along with their im- 
plications. 



5.1 Discovery of the extended bright structure in 
Orochi 

We find that Orochi is spatially resolved and that is has two 
components — an extended structure and a compact unre- 
solved one — from the analysis of the visibility amplitude 
as a function of the projected baseline length in SMA data 
(Fig. [S]). The plot of the phase calibrator is almost constant, 
whereas the visibility amplitude of Orochi appears to reduce 
significantly for a longer uv-distance. The visibility ampli- 
tude was fitted with a combination of a Gaussian (FWHM 
of 4".0) and an unresolved point source. This fitting was 
carried out using the uvfit task in MIRIAD. The estimated 
flux densities of the extended and compact structures are 
33.7 ± 16.2 mJy and 38.9 ± 3.2 mJy, respectively; in other 
words approximately a half of the total 880 /im flux (~73 
mJy) is originated from the extended component. Note that 
the flux densities of the extended component and the unre- 
solved component are 42.1 ± 20.3 mJy and 48.6 ± 4.0 mJy, 
respectively, after applying a correction of the primary beam 
responsE0. 

On the VLA image, Orochi is spatially resolved with 
a size of l".84 x 0".95 at PA 23.9°. A single source model 



^ In order to confirm the existence of the extended component, we 
carried out imaging of the extended component from the resid- 
ual data after subtracting the unresolved component in uv do- 
main. For this imaging, we used uv data at uv-distance = 10 — 30 
kA in order to achieve better signal-to-noise ratio. The image 
achieve the noise level of 4. 1 m Jy and the synthesized beam size 
of 4". 8x4". 3 (P.A.=84°). We found a point source on the residual 
image and its flux density is 20.8 ± 4.1 mJy. This flux density is 
consistent with the estimated flux density at uv-distance of 10-30 
kA on the visibility amplitude plot, 18.9 it 4.8 mJy (Fig. 6). In 
applying primary beam correction for Orochi, the effect of the 
extended source size should be considered, because Orochi is lo- 
cated 9 arcsec away from the phase centre of SMA. The primary 
beam response is 0.80 at the center of the source position, 0.87 at 
the inner edge of the extended components , and 0.73 at the outer 
edge of the extended components. This means that the effect of 
primary beam correction for the extended components depends 
on its structure. In order to evaluate this effect, we applied pri- 
mary beam correction for the residual image. Considering the fact 
that the extended component is unresolved in the beam size of 
this map (~ 4". 7), the corrected flux density on the residual im- 
age is equivalent to the corrected flux density of the extended 
component. The corrected flux density of the residual image is 25 
mJy, which is equivalent to applying the primary beam response 
of 0.8 , therefore we can apply the primary beam response of 0.8 
to the extended component. 
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Figure 5. The Z-Spec mm-wavelength spectrum of Orochi. The top panel shows observed flux density by Z-Spec. The dotted hne is the 
fitted continuum, S oc ;/3.00±0.li "pj^^ blue, red and green data points show the flux densities and error bars (without absolute calibration 
error) measured with the SMA, AzTEC/ASTE, and CARMA, respectively. The horizontal bars give the bandwidth. The yellow box 
shows areas used for calculating the continuum flux densities in Table 5. The bottom panel shows the Z-Spec spectrum after subtracting 
continuum flux. The yellow histogram shows the residual spectrum of Orochi and the blue histogram shows the 3 a errors. 



with this VLA source size can also fit the data, but the two- 
component model provides a better fit (Fig. [6|. The exact 
reason for the difference between the source size determined 
via our VLA and SMA data is no t clear. The VLA image 
published by ISimpson et all (|2006l ) was made with visibili- 
ties from the VLA's DnC and B configurations and yields 
a flux only 30 fiJy (~10 %, and < la) below the deeper 
VLA image used here (Arumugam et al., in prep), which 
supplements the Simpson et al. data with visibilities from A 
configuration. It is unlikely, therefore, that any loss of sig- 



nificant extended structure is due to the uv coverage. It is 
possible that the radio-to-submm fiux ratio varies across the 
structure: a dust temperature gradient may produce such a 
difference, because a lower dust temperature results in a 
lower radio flux density, based on the radio-FIR correlation 
(e.g., Condon 1992; Ivison ct al. 2010b) as shown in Fig. 12. 
It is also possible that faint, compact AGN-related emission 
may influence the radio morphology, leading to a smaller ap- 
parent size. It is clear that more sensitive submm and radio 
imaging is required to address the issue conclusively. 
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This two-component feature can also be observed in 
the CARMA 1300 ^Im image (Fig. [Tjl, although the signal- 
to-noise ratio of the CARMA data is worse than that of 
the SMA data. We conducted two-component fittings of the 
CARMA data as we did for the SMA data. The resultant 
source size of the extended component on the CARMA data 
is ~2".7 (FWHM). The estimated flux densities of the ex- 
tended and compact structures are 8.5 ± 3.5 mjy and 11.4 ± 
2.5 mJy, respectively, and the total flux density is ~20 mJy 
at 1300 /im (Table [T} . We also applied primary beam cor- 
rection for CARMA image using a primary beam response 
of 0.81, as we did for SMA image, because Orochi is located 
10 arcsec away from phase center on CARMA image. The 
corrected flux densities of the extended component and the 
unresolved component are 10.5 ± 4.3 mJy and 14.1 ± 3.1 
mJy, respectively. The estimated source size of the extended 
component by SMA and CARMA is also listed in Table 7. 
We suggest that the SMA and CARMA source sizes are 
consistent if we consider the uncertainties of the measured 
sizes. 

The discovery of the extended submm/mm bright com- 
ponent is distinguishing from almost all previously stud- 
ied SMGs. For instance, the median source size moderately 
bright SMGs (with a typical flux of ~2 mJy at 1300 or 
~7 mJy at 850 ^m) was measured to be ~0".4; no sourc e 
larger than l".2 has been reported (|Tacconi et al.ll2006l ). 
High-resolution radio observations are used to derive the 
median so urce size of 0".65 in S CUBA SMGs in the Lock- 
man Hole (|Biggs fc Ivisonll2008l ). Recently, brig hter SMGs, 
AzTE Cl (15 mJy in 890 /xm) and GN20 (27 mJv: llono et~al] 
l2006al ). have also been observed with SMA; their source sizes 
were estimated as app roximately 0".3 and 0".6, respectively 
l| Younger et al.l [20071 ): these values are similar to those of 
fainter SMGs. However, extended submm/mm structure at 
high redshift is observe d in HzRGs such as B3 J2330-I-392 7 
l|De Breuck et all |2003| ) and 4C 60.07 l|lvison et al.ll2008h . 
A lensed SMC, SMM J02399-0136, is known to have an ex- 
tende d mm bright structu re, and its corrected source size is 
~5". l|lvison et al.ll2010d ). 

It should be noted that uv data suggesting the presence 
of extended emission can also be fitted with a model with 
multiple sources. For instance, AzTECll, another bright 
SMG, can be model ed by two compact components with 
a separation of ~2" (|Younger et al.ll2009l ). However, in the 
case of Orochi, we found no multiple compact components 
in the VLA image with the beam size of l".3. Therefore, 
Orochi has certainly an extended component, although the 
size has not been precisely determined yet. 



5.2 SED and redshift of Orochi 

We conducted SED fitting and estimated the photometric 
redshift of Orochi using two methods. One involves the fit- 
ting of an SED to the optical and NIR data, and the other 
involves the use of submm/mm and radio data. The former 
corresponds to fiux from stellar emission, and the latter, to 
cold dust emission. 



"18'"31?2 31 TO 



30r8 



30f6 



30f4 



30r2 



Right Ascension (J2000) 

Figure 8. The SMA 880 ^tm image of the extended component of 
Orochi, produced by the residual visibility data after subtracting 
unresolved component in the uv domain (See section 5.1). The 
contour levels are 2.5, and 5 cr {I a = 4.1 mJy). The synthesized 
beam is shown on the lower left. The white triangle marks the 
position of the compact component. 



Table 7. Summary of submillimeter /millimeter bright compo- 
nents. 



Size 



Flux density 



SMA 880 



compact unresolved 48.6 ± 4.0 mJy 
extended 4".0 ± 1".2 42.1 ± 20.3 mJy 



CARMA 1300 ^lm 



compact unresolved 
extended 2". 7 ± 1".2 



14.1 ± 3.1 mJy 
10.5 ± 4.3 mJy 



5.2.1 Optical and NIR photometric redshift 

The optical/NIR photometric redshift is calculated using the 
code Hyperz (|Bolzonella et al]|2008l ). with ur, B, V, Rc, i , 
z , J, and K. The default SED template of IColeman et al.l 
l|l980l ), extended using the lBruzual fc CharlotI (|l993l ) model 
is adopted. The obtained redshift is 1.39 (x^ ~ 4.7), with 
the 68% confidence interval of 1.383 - 1.402. 

Next, by fixing the redshift to 2 = 1.39, we 
carried out SED fitting using the GALAXEV model 



(|Bruzual fc CharlotI 12003 ) to obtain its physical properties 
of the galaxy. Two star formation histories of instantaneous 
burst and continuous burst models a re adopted, wit h the 
Salpeter initial mass function and the ICalzettil (|200G| ) dust 
extinction curve. The results are shown in Fig. |9] it can 
be observed that the optical/NIR SED is well fitted by the 
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Figure 6. 880 visibility amplitude as a function of the baseline lengths (uv distances) toward Orochi (left) and the visibility calibrator 
(right), along with the results of Gaussian and point source fitting. These plots only show the real part of the visibilities, and the flux 
densities are not primary-beam corrected on these plots. The visibility amplitude of the calibrator (right) is well modeled with an 
unresolved point source, whereas Orochi can be fitted with two components, i.e., a combination of a Gaussian shaped component with a 
FWHM of ~4" (dash-dot-dot-doted line) and an unresolved point source (dotted line). The resultant fitting is indicated by a dot-dashed 
line. The solid curve line shows the amplitude function derived from the VLA source size (1".84 X 0".95 at PA 23.9°), and its amplitude 
is scaled to be fitted to the SMA result. 
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Figure 7. 1300 /.tm visibility amplitude as a function of the baseline lengths (uv distances) toward Orochi (left) and the visibility 
calibrator (right), along with the results of Gaussian and point source fitting. These plots only show the real part of the visibilities, and 
the flux densities are not primary-beam corrected on these plots. The visibility amplitude of the calibrator (right) is well modeled with 
an unresolved point source, whereas Orochi can be fitted with two components, i.e., a combination of a Gaussian shaped component 
with a FWHM of ~2".7 (dash-dot-dot-dotcd line) and an unresolved point source (dotted line). The resultant fitting is indicated by a 
dot-dashed line. 



500 Myr instantaneous burst model, with a stellar mass of 
1.8 X 1O"M0 and a dust extinction of E{B -V)= 0.31. 

5.2.2 submm/mm and radio photometric redshift 

We used a Monte Carlo photometric redshift method 
ijHughes et all I2OO2I : lAretxaga et all l2003l ) to estimate the 
phototmetric redshift based on submm/mm and radio data. 

A Monte Carlo photometric redshift method has been 
developed to take into account constraining prior infor- 
mation such as a range of local SED templates, the 



favoured evolving luminosity function of dust enshrouded 
star-forming galaxies to 2 « 2, and the amplification of cer- 
tain fields. We assume that the SEDs of SMGs are well repre- 
sented by 20 SEDs of local starbursts, ULIRGs, and ACN to 
provide FIR-radio SEDs. These SEDs cover a wide range of 
FIR luminosities (10 ^ ° - 10^^-^ Lq) and te mperatures (25- 
65 K). See details in lAretxaga et al] (|2007l ) for the method. 
We adopt a very conservative estimate for the amplification 
of the source (/i = 10). 

We used the FIR to radio portion of the SED, i.e., the 
photometric data at 160 /im (Spitzer; upper limit), 1000 
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Figure 9. Result of Optical/NIR photometric redshift and SED fitting of Orochi. Left panel shows the reduced versus photomet- 
ric redshift obtained by Hyperz code. Right panel shows the best-fit optical/NIR synthetic model spectra, by fixing the redshift to 
1.39 obtained by Hyperz. The dotted and dashed lines show the resultant spectra assuming instantaneous and continuous starburst, 
respectively. 



fim (Z-Spec), 1100 fim (AzTEC/ASTE), 1500 /im (Z-Spec) 
, 20 cm (VLA) and 50 cm (GMRT), for the derivation of 
the photometric redshift. These three photometric data at 
submm/mm wavelength were selected as representative s of 
the Rayleigh- Jeans slope. 

The derived redshift probability distribution is shown 
in Fig. 9. The estimated redshift is SA'^q I^. The peak of the 
probability density distribution is very stable to changes in 
magnification. Solution were derived FIR amplification fac- 
tors ranging from 1 (no magnification) to 50, finding peak 
values between 3.3 and 3.5. As the magnification increases, 
we find a slight transfer of the probability density towards 
higher redshift values, which displaces the 68 percent confi- 
dence interval up to 3.41q;5. 

5.2.3 z ^1.4 or ^3.4 ? 

We now discuss the inconsistency between the two photo- 
metric redshift results. The optical/NIR photometric red- 
shift suggests z ~ 1.4, while the submm/mm and radio pho- 
tometric redshift suggests z ~ 3.4. Next, we show that this 
apparent discrepancy can be understood by observing an 
optically dark SMG lying at z ~ 3.4 behind a foreground 
galaxy located around z ~ 1.4, visible at optical/NIR wave- 
lengths. 

The optical/NIR-based photometric redshift of ~ 1.4 
seems to be robust because we obtain nearly continuous 
detections in optical/NIR, i.e., ~0.4 /im-2.2 /im; its SED 
shows a clear break feature at around 1.2 ^m, which is 
very likely to be a 4000 A break around 2 ~ 1.4. Although 
the IRAC colours of Orochi, log(S5.8fim/5'3.6,jm) = 0.21 and 
log(5'8.0mn/»5'4 Sum) = 0-28, su ggest the existence of an ob- 
scured AGN (|Lacv et al.ll2004h . such a sharp 4000 A break 
can be reproduced solely by stellar spectra, (i.e., it is very 
difficult to produce a clear break with AGN models) . There- 
fore, we suggest that the observed optical/NIR spectrum 
is dominated by a stellar component; Further, the possible 
4000 A break feature reliably constrains the redshift. The 



resultant reduced distribution shows a sharp peak at z 
~ 1.4. 

On the other hand, the submm/mm and radio photo- 
metric redshift is also fairly well constrained with almost 
continuous coverage from 880 /im to 1500 /im combined with 
radio and MIPS upper limits. Although it is possible to fit 
the observed SED with an SMG at z — 1.4, it requires an 
SED template with a very low dust temperature, Tdust ~ 
20 K, which is nearly the lowest measured dust temp erature 
distribution among SMGs (e.g., H ovdcs et al.l 12009 ). Refer 
to the next section for a detailed discussion of SED modeling 
with a grey-body. 

Another constraint on the redshift of the gas/dust com- 
ponent comes from Z-Spec observations. In general, the line- 
to-continuum ration (L/C ratios) are expected to be smaller 
for higher J transitions because the continuum flux increases 
with On the other hand, the CO line fiux increases 

with or for thermalized conditions, and since high-J 
lines are often subthermally excited, they are much weaker 
than in a thermalized case. 

As no emission line features are detected for the intense 
continuum emission, the upper limit on the L/C ratio is 0.1- 
0.3 (2 a) accross the Z-Spec band. This value is low if we 
assume a redshift of ~ 1.4, where CO(J=4-3), CO(J=5-4), 
and CO(J=6-5) whould be expected to be in the Z-Spec 
band. 

We collected the L/C ratios from the measured CO 
fiuxes and adjacent continuum emission using interferom- 
eters in the literature and we found that L/C ratios are typ- 
ically in the range of 3-20 for high-z quasar s (pDownes et al. 
1999'; 'WeiB et al' '2007'; ' Omont et al. 1 L1996'; Bradfor d et al.l 
2009; Wang et al. 2010 ), radio galaxies (J-'apadqpoulqs et al. 
2OO0i) . SMGs ifCreve et al]|2005l : iTacconi et al.1 120061 1^1. and 



^ There are no direct descriptions of the continuum fluxes, and 
L/C ratios are estimated from the CO(J=4-3) spectrum in these 
papers. 
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Figure 10. The left panel shows the colour-colour-redshif plot for Orochi. The flux ratios of the mock galaxies are represented as 
diamonds, and their redshifts are colour-coded according to the scale shown in the panel on the right. The cross represents the measured 
colours of Orochi, and the dashed box shows the 1 a uncertanity in each colour. The top right panel shows estimated redshift probability 
distribution of Orochi. In the bottom right panel, the observed SED of Orochi normalised to the flux densities at 1000 lira is shown as 
squares and arrows. The arrow indicates 3 a upper limits. The squares denote the detection at a level ^ 3 a, with 1 a error bars. The 
template SEDs (lines) are redshifted to z = 3.4. The template SEDs at this redshift compatible within 3 a error bars with the SED of 
Orochi are displayed as blue lines. The photometric data are shown at 1000, 1100, 1500 fim, 20 and 50 cm, and The upper limit is shown 
at 160 fim. 



a local ULIRG ([Matsushita et aLllioogh . This indicates that 
we have already achieved sufficient line sensitivity for the 
observed continuum fiux level if Orochi is lying at around 
z ~ 1.4. 

On the other hand, the observed low L/C ratio is ac- 
ceptable if Orochi is at 2; ~ 3.4; this is because the observed 
J lines such as CO(J=9-8) to CO(J=13-12) in the Z-Spec 
band are considered to be of a much higher order. Smaller 
L/C ratios around unity h ave been reported f o r CO(J=9- 



in APM J08279-h5255 llDownes et al.l 



[2003) and the Cloverleaf (jBradford et al 



999I: IWeifi et~all 



2009t ). It should 



be emphasized that these two quasars show e xceptionally 
highl y excited CO lines up to CO(J=10-9) level (jWeifi et al.1 
the measured CO SED suggests that the CO fluxes 
are consistent with the thermalized condition as they are ex- 
pected to increase with up to J ~ 10. Other typical high-z 
quasars as well as starburst cores of local IR bright galaxies 



like M82 and NGC 253 become subthermal at around J ~ 
4 - 6; a simple LVG model calculation suggests that the ex- 
pected C0( J=9-8) flux in a typical starburst core (e.g., gas 
density ~ 10* cm"'' and kinetic temperature Tkin = 60 
K) will be smaller than the thermalized-case flux by a fac- 
tor of 20 or more; i.e., it is natural to expect that the L/C 
ratio for these high- J lines are much smaller than unity, as 
observed in Orochi with Z-Spec. 

Here, we briefly comment on the constraint on redshift 
from [CII] emission (frcst = 1.9019 THz); this line can be 
in the Z-Spec band if 5.2 < z < 9.5. The upper limit to 
the [CII] line luminosity is Lfcn] < 3.6 x W^Lq (2 a), if 
it is at z = 5.5 with a line widt h of 780 km s~^, which 
is a median value for SMGs (e.g., ICreve et al.ll2005l ). This 
results in a strict upper limit on the [CII] to FIR lumi- 
nosity ratio, L[cii]/iFiR < 3.6 x 10~® (2 a), which is al- 
ready much smaller than the observed L[cii]/^fir ratios 
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for high-z quasars/SMGs, i .e., ^10 — 10 ^ (jlono et al.l 
l2006bl : iMaiolino et al.ll2009l : llvison et al.ll2010al '). We, there- 
fore, suggest that the redshift of Orochi would be less than 
~ 5. 

In summary, we suggest that the mm/submm bright 
component of Orochi is likely to be at 2 ~ 3.4, whereas 
the optical/NIR counterpart candidate with a photometric 
redshift of ~ 1.4 is a foreground galaxy with an old stellar 
population. Although it is difficult to ignore the possibility 
that both systems, the i.e., submm-to-radio bright compo- 
nent and its optical/NIR counterpart, are lying at the same 
redshift (i.e., at 2; ~ 1.4), this possibility requires an excep- 
tionally low dust temperature for SMGs (Tdust ~ 20 K) as 
well as a low-excitation molecular gas for which the mid-J 
lines such as CO(J=4-3) are not well excited. 

In the following sections, we estimate physical quanti- 
ties for two cases — 2 ~ 3.4 and 2 ~ 1.4 — because it 
is difficult to completely reject the latter possibility at this 
moment, even though the former seems more probable. 



5.2.4 Dust SED model using a grey-body 

We model SEDs from FIR to radio wavelengths using grey 
bodies. To draw SEDs at radi o wavelengths , we assume 
the ra dio-FIR correlation (e.g.. [CondonI [l992l : llvison et al.l 
l2010bl ). This correlation is described as the q- value Q It is 
known t hat the mean of qtir for SMGs is 2.32 with a scatter 
of 0.34 (jMichalowski et al.llioiol '). We use this typical qtir 
of known SMGs for extrapolating the radio spectrum. More- 
over, LiR is derived from the SED of a model grey-body. 

We also derived the radio emissivity index a [S (x 1'°') 
of Orochi. Orochi has two photometric data at radio wave- 
lengths, 20 cm (VLA) and 50 cm (GMRT). From these 
data, the radio emissivity index a of Orochi is derived as 



p=1 .4, T=40 K 
P=1.5, T=30 K 
p=1 .9, T=20 K 
SXDF1 100.001 



+0.29 



-0.33 
0.7 to -0.8: 



Comparing a of local star-forming galaxies 
Condon. 1992, ), and SMGs at 2 ~ 2 (-0.75 ± 0.06; 



llbar et aLllioiOl ). the measured a in Orochi is found to be 
rather closer to a pure thermal radio spectrum. If Orochi is 
located at 2 > 3, the a of Orochi seems to be consistent with 
predictions that non-thermal radio continuum emission will 
be suppressed at high 2 (2 > 3 or so) due to the increased 
energy losses from the inverse Compton scattering off in the 
cosmic microw ave background to the cosmic ray electrons 
l|Murphvll2009l ). 

For the models based on the optical/NIR photometric 
redshift of 1.4, the dust emissivities were derived as /3 = 1.4, 
1.5, and 1.9 for Tdust = 40, 30, and 20 K, respectively, by fit- 
ting the Z-Spec continuum data at 1000-1500 pim. We find 
that a low dust temperature, 20 K, is required if a photo- 
metric redshift of 1.4 is assumed. Possibilities of higher dust 

^ The q-value is defined as g oc log(LpiR) — log(7i.4GHz), where 
7i 4QHz is a rest-frame 1.4 GHz luminosity density and Lpm 
is defined as the 42.5-122.5 /^m FIR luminosity. This relation 
corresponds to the median q-value, < q > ~ 2.3, with a rms 
scatter of ^0.2 in nearby normal gal axies dC ondon ' 1 9921 ) and 
2.34 ± 0.01 in local IRAS galaxies l lYun et a l. 2001,). On the 
other hand, we use the total 8-1000 /xm IR (Lir) as an indi- 
cator of dust emission. In this case, the q-value is redefined as 
ijTiR <x log(LiR) - log(/i,4GHz). The value of ijtir is 2.64 ± 0.02 
in loc al galaxies, which is consistent with the above q-value llBelll 
l2003t) . 




Figure 11. The SED of Orochi and SED fitting results in FIR- 
radio wavelengths, assuming z=l.A. We calculated grey-bodies 
for T^uat = 20, 30, and 40 K. The /9 for each grey-body is es- 
timated by fitting the Z-spec continuum data at 1000-1500 /^m. 
The dashed lines show the scatter of the q-value. The upper limits 
are shown at 24, 70, and 160 lira, and the photometric data are 
shown at 880, 1000, 1100, 1200, 1300, 1400, 1500 /xm, 20 and 50 
cm. We find that the dust temperature of Orochi must be as low 
as 20 K if we assume z ~ 1.4. 



temperatures such as 30 K and 40 K are clearly eliminated 
because the upper limits at 24, 70, and 160 /xm bands tight 
constraints on the dust temperature. The radio flux density 
also seems to be consistent with a model with Tdust of 20 
K. Again, higher dust temperatures such as 30 K and 40 K 
are not acceptable because at such high dust temperatures, 
LiR is overluminous in terms of radio luminosity, as clearly 
shown in Fig. 1111 

On the other hand, for models based on the 
submm/mm-to-radio photometric redshift of ~3.4, we find 
that the observed submm/mm to radio SED can be well ex- 
plained by the dust temperature of ^ 30 K, as shown in Fig. 
12. This is consistent with the typical dust temperature of 
SMGs (e.g.. lKovdcs et al.ll2006D . 



5.3 A possible positional offset between a 
submm-radio bright component and an 
optical/NIR counterpart candidate 

Here, we present a detailed comparison of the multi- 
wavelength images of Orochi. Fig. 13 displays the peak po- 
sitions of CARMA, SMA, and VLA sources listed in Ta- 
ble 3 superposed on the Subaru 2'-band image. We find 
a possible positional offset between the optical peak and 
the submm/mm/radio peaks; its separation is approxime- 
tely 0".l-0".7. Although these separations are smaller than 
the beam sizes of the observations, the evaluated combined- 
astrometry errors listed in Table 3 are smaller than the 
claimed offset, especially in declination. Furthermore, all 
three observations (i.e., CARMA, SMA, and VLA) show 
a similar tendency, i.e., the peaks of these wavelengths 
are located at the north west part with respect to the 
2'-band peak, which is also associated with the K band 
peak. These facts suggest the possible existence of a spa- 
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P=1.6, T=50k 
P=1.8, T=40 K 
P=2.0, T=30 K 
Orochi 




10'' 10 10= 

Figure 12. The SED of Orochi and SED fitting results in FIR- 
radio wavelengths, assuming 2=3.4. We calculated grey-bodies 
for T^ust = 30, 40, and 50 K. The 13 for each grey-body is es- 
timated by fitting the Z-spec continuum data at 1000-1500 /im. 
The dashed lines show the scatter of the q- value. The upper limits 
are shown at 24, 70, and 160 fim, and the photometric data are 
shown at 880, 1000, 1100, 1200, 1300, 1400, 1500 /xm, 20 and 50 
cm. We find that the dust temperature of Orochi seems to be as 
low as 30 K if we assume z ~ 3.4. 



' I ■ ■ ■ I ■ ■ ■ I ■ 



30.74 2:18:30.70 30.66 
Right Ascension (J2000) 



I 



24 


u 

0) 




(0 




n 


22 






< 




Dl 




ta 


20 


S 


18 





Figure 13. Radio/mm/submm peak positions superposed on the 
z-image of Orochi. The white, black, and blue crosses (express- 
ing the position uncertainties) correspond to the SMA, CARMA, 
VLA peak positions, respectively. The unit of the background 
z'-band image is Mag (AB) per arcsec^. 



tial oflFset between the peaks in the optical/NIR and the 
submm/mm/radio wavelengths. 

As discussed in section 5.2.3, the optical/NIR bright 
counterpart candidate can be a foreground source of Orochi. 
In this case, the extended bright structure observed in the 
optical and NIR bands, i.e., the north-west extension from 
the optical/NIR peak positions around the SMA, CARMA, 
and VLA is probably a true optical counterpart of Orochi 
and is blended with the foreground source. 

In order to confirm this possibility, we need higher an- 
gular resolution images for all wavelengths. Future ALMA 
observations combined with HST or ground-based AO ob- 
servations will be conducted for addressing this issue. 



5.4 The Origin of 1100 Flux in Orochi 

The observed 880 fim flux of Orochi, ~ 91 mjy, makes 
this one of the brightest SMGs known to date. For in- 
stance, it is much brighter than the bright est SMGs re- 
veale d by SHADES (~ 10-15 mJy at 850 /xm; ICoppin et al] 
I2OO6I ) and brighter or comparable to the known lensed 
submillimeter bright high-z sources such as the Clover- 
leaf (58.8 ± 8.1 mJy at 850 fim , with a lens magnifi - 
cation of 11; iBarvainis et all |2002| : IVenturini et al.ll2003l ). 
APM08279-I-5 225 (84 ± 3 mJy at 850 /xm, with a lens mag- 
nification of 7; IBarvainis et al.ll2002l : iLewis et al.ll2003 ) and 
SMMJ2135-0102 (1 06.0 ± 7.0 mJy at 87 /im, with a lens 
magnification of 32: ISwinbank et al.ll2O10l ). A similar bright 
SMC has also been uncoverd by MA MBO on IRAM 30 m 
observations (30 ± 2 mJy at 1200 fim; iLestrade et"al]|2010D . 
SPT surveys detected 47 dust-dominated mm sources above 
10 mJy at 1.4 mm in 87 de g^ , and 20 of 47 do not have coun- 
terparts at low-z universe (|Vieira et al.ll20l"ol ). If all of 20 are 
indeed SMCs located at high-z universe, a surface density of 



these ultra-bright SMGs is ~ 0.23 sources per 1 deg^. Con- 
sidering we found only 1 ultra-bright SMG in about 0.33 
deg^ , which is surveyed area in SXDF, it is consistent that 
Orochi belongs to the same population with SPT sources. 

It is possible that Orochi is gravitationally lensed, be- 
cause many of the ultra-bright SMGs are often strongly 
amplified by a foreground massiv e cluster of galaxie s (e.g., 
Wilson et al. 2008b; Johansson o t al.ll2010l : iSwinban k et al.l 
I2OI0I : llvison et al.l l2010al ). The recent model analysis for 
submm/mm count shows that these ultra bright SMGs are 
likely to be gravitationally lensed by foregrou nd massive 
halos (|Lima et al.ll2010l : iHezaveh fc Holdeij|2010l '). Ahhough 
we found no catalogued foreground clusters around Orochi, 
analysis of photometric redshifts of optically visible objects 
around Orochi suggest that there are cluster candidates at 
z ~ 0.4 and 1.4, and they may be partially responsible for 
lensing. 

It is also possible that this SMG is gravitationally lensed 
by other foreground galaxies near Orochi. Theoretically, 
when a galaxy causes a lensing effect on another galaxy, the 
former galaxy is located at around less than 1 arcsec from 
the latter. In the case of Orochi, the optical/NIR bright fore- 
ground galaxy at z ~ 1.4 can be a candidate for the lensing 
galaxy. The estimated stellar mass of the optical/NIR bright 
source by SED fitting is ^ 10^^ Mq; therefore, it is suffi- 
ciently massive to cause gravitational lensing. In fact, some 
of 500/im-selected ultra bright SM Gs are indeed gravita tion- 
ally lensed by foreground galaxies (jNegrello et al.llioTOi ') . 

Nevertheless, we note that there is no clear evidence 
for strong lensing in existing high-resolution optical/infrared 
images, for instance, the presence of multiple and/or dis- 
torted structures. This situation is different from that of 
SMM J2135-0102, another uhra-bright SMG, where dis- 
torted/elongated images can be clearly observed even in the 
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resolution of the IRAC bands, although the noise level of 
the IRAC bands in SXD F is at least 10 times better than 
that of SMM J2135-0102 (|Swinbank et al.ll201(]| ). 

On the other hand, the situation of Orochi is similar to 
that of HDF850.1 (7.0 ± 0.4 mjy at 850 fj.m- Hughes ct al. 
1 19981). The photom etric redshift of HDF850.1 is z = 4.1^0.1 
l|Cowie et al] l2009t ) and it is indicated that HDF850.1 is 
gravitationally lensed by an elliptical galaxy at z =1.22 
about 1 arcsec away from HDF850.1 (|Hughes et all 1 19981 : 
iDunlop et al. 200l: Icowie et al.ll2009h . HDF850.1 also has 
no clear counterpart image in the mid- and near-infrared 
deep images. 

The elevated mm/submm fluxes can be from non- 
thermal synchrotron emission if they are powered by a radio- 
loud AGN. In fact, wide-fiel d mm surveys often uncover 
such non-thermal sources (e.g.. lVoss et al l2006l : IVieira et al.1 
I2OIOI ). In the case of Orochi, however, we conclude that the 
submillimeter flux is dominated by thermal dust emission 
heated by massive stars and not by a radio-loud AGN be- 
cause the submm-to-radio SED can be well reproduced by 
the known LiR-to-radio correlation observed in previously 
studied SMGs and local starburst galaxies, as shown in sec- 
tion 5.2.4. Absence of the time variability in the submm/mm 
fluxes of Orochi, i.e., in the 1100 pm flux measured by 
AzTEC in November/December 2008, the 1300 ^lm flux by 
CARMA in August 2009, the 880 ^m flux by SMA in De- 
cember 2009, and the Z-Spec data in November 2009, may 
also suggest the thermal origin of the elevated flux of Orochi. 

5.5 Estimation of physical properties of Orochi 

We estimate the physical quantities of Orochi, i.e., Lir, dust 
mass (Md) and gas mass (Mgas), by adopting two possi- 
ble redshifts and models in section 5.2: (1) z = 1.4, for 
which the SED model of the grey-body with = 20 K 
and P — 1.9, and (2) z = 3.4, for which we used an av- 
eraged SED template of the 4 best fltted SEDs (the SED 
templates of NGC3227, NGC7771, IRAS05189-2524, and 
IRAS12112-h0305), providing =35 K and /3 = 2.0. Here, 
we simply show the physical quantities of Orochi using the 
apparent luminosity, i.e., Lir = M^iri where Ljr is the in- 
trinsic luminosity and ^ is the magniflcation factor, because 
we have no constraint on the magnification factor (if any) 
at this moment. Furthermore, here, we do not distinguish 
between the compact and extended components. 
We estimated SFR using the equation below. 

SFR{MQyi'^) = 4.5 x 10~''''LiR(ergs s~^) (2) 

ljKennicutt.lll998t ). Md is derived as 

(1 + 2) Kd[l'rest)B{Vrest,Tdust) ' 

where Soba is the observed flux density. In this paper, the 
flux density is at 1100 /im, i^rest is the rest-frame fre- 
quency; Kd{vrest), the dust mass absor ption coefficient; an d 
B{i'rest,Tdust), the Planck function (iHughes et al.l [l997l ). 
We assume that the absorption coefficient varies as Hd . 
We adopt Kd(125 urn) = 2.64 ± 0.29 m^kg ~\ the average 
value of various studies (|Dunne et ahllioOa ). The gas mass 
is then derived using a gas-to-d ust mass ratio of 54 , which 
is an average value for SMGs in iKovdcs et al.1 (120061 ). 

The derived quantities are listed in Table [8] along with 



the gas consumption time scale, which is defined as the ratio 
of Mgas/SFR or 1/SEF, where SEE refers to the star forma- 
tion efficiency. However it should be noted that many of 
these results are rough estimates based on median relations, 
assumption about the SED. 

Table |8] shows that if 2 ~ 3.4, Orochi will consume its 
gas reservoir within a short time scale (~ 3 x lO'^ yr), which 
is indeed comparable to thos e in extreme star bursts like the 
hearts of local ULIRGs (e.g.. iKennicutt .11 19981 ) . On the other 
hand, if this source has a lower 2 of ~ 1.4, Orochi has a fairly 
long gas consumption time scale of '^l xlO^ yr, which is 
almost comparable to that of star- forming clouds in the disks 
of spiral galaxies in the local universe, and hence, rather 
unlikely. 

5.6 Starburst form of Orochi 

One of the remaining key questions of Orochi is its form 
of starbursts: What type of star formation or starburst pro- 
duces such an SFR elevated to ~ 870 Mq yr~^ (if 2 ~ 1.4) or 
11000 Mq yr"'^ (if 2 ~ 3.4)? Recent high angular resolution 
studies of high-z quasars and SMGs suggest that Eddington- 
limited maximal starbursts often occur at the heart of thes e 
extreme objects (lYounger et al.|[2"008l : iRiechers et al.ll2009l ). 
In these galaxies, the peak SFR and Lir surface densities 
reach up to Espr ~ 10^ Mq yr"^ kpc"^and Eljj^ lO" 
Lq kpc~'^, respectively; these values are close to the the- 
oretical maximum value imposed by the Edd ington limit 
l|Elmergreen et al.lll999l : [Thompson et al.ll2005l ). Therefore, 
it is intriguing to address whether the SFR/Lfir surface 
densities at the center of Orochi are close to the maximum 
starburst condition or not. 

We estimate the surface gas mass density (Sgas), SFR 
surface density {Y,sfr) and Lir surface density (S^jp^) for 
these two components (Table 7), respectively. 

We assume that the source size of the unresolved com- 
pact component is 0".4, which is t he median size of the 
submillimeter brigh t region of SMGs l|Chapman et al. I l2004l : 
iTacconi et al.|[2006l ). The derived surface densities for two 
possible photometric redshifts, ~ 1.4 and --^ 3.4, are listed 
in Table H 

The Tigas and other derived densities of the extend com- 
ponent are similar to or higher than those o f the star-forming 
disk regions in local spiral gal axies (e.g., iNishivama et al] 
I2OOII : I Wong fc Blitz et al.|[2003 ). 

This indicates that the extended component of Orochi 
is a large lump of relatively moderate starburst regions, al- 
though we should note that the extended structures can be 
a lensed image. We then need to estimate the magnification 
factor to discuss the quantitative properties of star forma- 
tion. 

On the other hand, the derived parameters of the com- 
pact component are rather close to those of ULIRGs in 
the local universe, i .e., the predicted maximum limit (e.g., 
IRiechers et al.l [20091 ) . although it still depends on the as- 
sumed source size and magnification factor (if any). 



SUMMARY 

An ultra-bright SMG, Orochi, has been detected using with 
the AzTEC on ASTE. Subsequent CARMA, SMA, and Z- 
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Table 8. Estimated IR luminosity, dust mass, gas mass, star formation rate, and inverse of star formation efficiency (corresponding to a 
gas consumption time scale) for two possible redshifts, z = 1.4 (model 1) and 3.4 (model 2), which are esitimated in section 5.2.4. Here, 
we describe the observed quantities as M^jj^, M^^^, fiM*^, and /iSFR* , where indicates the magnification factor. Note that 1/SFE is 
independent of magnification. 



model z p ^iL*^ fiM* M^^gas fiSFR* 1/SFE 

(1) 1.4 20 K 1.9 5.1 X 10^2 I.TxIO^'^Mq 9.1 x lO" Mq 870 Mgyr"! 1.0 X lO'^ yr 

(2) 3.4 35 K 2.0 6.3 x lO^^ Lq 5.8 x 10^ Mq 3.1 x lO^i Mq 11000 Mgyr"! 2.8 x lO'' yr 



Table 9. Estimated surface density of molecular gas mass, SFR, and Li^ ( Sgas, Ssfr, and Sljj^, respectively) for two components 
revealed by SMA, i.e., spatially unresolved, compact component (the size was assumed to be 0".4) and spacially extended component 
(4".0). We adopted a scale of 8.432 kpc/arcsec for z = 1.4 and 7.396 kpc/arcsec for z = 3.4. Adopted SED models, dust temperatures, 
and parameters are the same as Table |8] 



Model 


z 


Tdust 


Component 


Size 


2, 




SsFR 








(1) 


1.4 


20 K 


compact 
extended 


0".4 (3.4 kpc) 
4".0 (33.7 kpc) 


5.4 x lO^o 
4.1 X 10* 


M0kpc~2 
M0kpc-2 


5.2 xlO^ M0yr~ 
3.9 xlO-l M0yr- 


ikpc-2 
-lkpc-2 


3.0 X 10" L0kpc^ 
2.3 X 10^ L0kpc- 


-2 

-2 


(2) 


3.4 


35 K 


compact 
extended 


0".4 (3.0 kpc) 
4".0 (29.6 kpc) 


2.4 X 10^0 
6.7 X lO'^ 


M0kpc-2 
Af0kpc-2 


8.3 xlO^ M0yr- 

7.4 xlO" M0yr- 


lkpc-2 
lkpc-2 


4.8 X 10^2 L0kpc- 
4.3 X lO^'' L0kpc- 


-2 
-2 



Spec observations confirm the AzTEC/ASTE detection of 
Orochi. The major findings and conclusions are summarized 
Eis follows. 

(i) We discovered a 37.27 ± 0.65 mjy source at A = 1100 
^m in SXDF using AzTEC mounted on ASTE. 

(ii) CARMA A =1300 fim and SMA A =880 /xm obser- 
vations successfully confirmed the AzTEC/ASTE detection 
of Orochi. The peak positions of the CARMA/SMA sources 
coincide with the AzTEC/ASTE peak position. In addition, 
these flux densities from 1300 to 880 ^m are consistent with 
a single SED for ^ = 1 ~ 2. The 880 fim flux density of 
Orochi is ~91 mJy, which makes it one of the brightest 
SMGs after SMMJ2135-0102 (~106 mJy at 870 urn). 

(iii) The CSO 10 m telescope equipped with Z-Spec 
was used to conduct a blind search for redshifted molec- 
ular/atomic lines and continuum emission in the 190-308 
GHz band. The measured continuum flux is consistent to 
the SMA/AzTEC/CARMA measurements, although no sig- 
nificant emission/absorption line features were found. The 
derived upper limit to the line-to-continuum flux ratio (L/C 
ratio) was 0.1-0.3 (2 a) across the Z-Spec band. 

(iv) We flnd that Orochi is spatially resolved and has 
two components, i.e., an extended structure (FWHM of 
~ 4") and a compact unresolved one, based on the anal- 
ysis of the visibility amplitude as a function of the pro- 
jected baseline length in both SMA and CARMA data. Ap- 
proximately half of the total millimeter/submillimeter flux 
arises from the extended component. The discovery of the 
extended submm/mm bright component is distinguishing 
because such an extended bright structure has not been ob- 
served in previously studied normal SMGs, which show a 
median source size of 0".4. 

(v) Multi-wavelength counterparts of Orochi have been 
identified near the peak position of CARMA/SMA using 



optical (Subaru), NIR (UKIRT; Spitzer), MIR(Spitzer) and 
radio (VLA and GMRT) images. 

(vi) A robust photometric redshift of 1.4 was derived us- 
ing optical/NIR data and the code Hyperz based on a detec- 
tion of possible 4000 A break feature. On the other hand, a 
photometric redshift using submm/mm and radio data sug- 
gests 2 ~ 3.4, eliminating the possibility of z < 3. A fairly 
low dust temperature, Tdust ~ 20 K, is required if z — 1.4, 
whereas Tdust of 30-50 K, typical for SMGs, can explain the 
observed FIR to radio SED for a redshift range of 3-5. 

(vii) The discrepancy in photometric redshifts can be un- 
derstood for an optically dark SMG lying at 2 ~ 3.4 with a 
foreground galaxy around z ~ 1.4. Indeed, we flnd a posi- 
tional offset of ~0".2-0".7 between sub/mm/radio peaks and 
and optical peak, implying a possible coincidental overlap of 
two objects along the line of sight. 

(viii) If Orochi is located at 2 ~ 3.4, higher order J-lines 
such as CO( J = 9 - 8) to CO(J = 13 - 12) faU into Z-Spec 
band. However we did not achieve any significant detection 
of lines. The upper limit on L/C ~ 0.1-0.3 in Z-Spec band 
is expected for these very high J lines. 

(ix) If the millimeter/submillimeter bright component of 
Orochi is indeed lying at 2 ^ 3.4, the deduced apparent 
FIR luminosity (Lir) and star formation rate (SFR) are 
LiR - 6 X 10^^ Lq and SFR ~ 11000 M0yr"\ respectively, 
if the huge Lir is originated from a massive starburst. The 
apparent surface densities of Lir and SFR, Egas and ^Lj^, 
of the unresolved compact component are similar to those of 
local ULIRGs cores, i.e., close to the theoretically expected 
maxmum value imposed by the Eddington limit, although a 
constraint is required on the magnification factor (if any) to 
yield further quantitative discussions. 
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